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Abstract: We construct contact interactions for bosonic and fermionic point par¬ 
ticles. We first relate the resulting theories to classical electrostatics by taking func¬ 
tional averages over worldlines whose endpoints are hxed to charged particles. Count¬ 
ing those paths which pass through a space-time point gives the static electric 
held at that point, provided we take the limit where the length measured along the 
worldlines is large. We also investigate corrections to the classical held that arise 
beyond leading order in this limit before constructing a theory of point particles 
that interact when their worldlines intersect. We quantise this theory and show that 
the partition function contains propagator couplings between the endpoints of the 
particles before discussing how this is related to the worldline formalism of quantum 
held theory and general action at a distance theories. 
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1. Introduction 

Classical electromagnetism is conventionally described by Maxwell’s field theory and 
there seems to be little room for debate about its formulation. In [1] and [2], however, 
building upon [3] an alternative approach to determining the electromagnetic field 
strength tensor for a pair of charged particles led directly to a novel interacting 
string theory. This theory contained contact interactions on the string worldsheet 
which served to produce expectation values of Wilson lines in Abelian quantum field 
theory. In the case of electrostatics, the description given in [3] was in terms of point 
particles whose worldlines have their endpoints fixed to the charged particles. The 
electric field at a position in space-time was arrived at via a weighted average over all 
such worldlines which also pass through the given point. The physical picture which 
motivated this approach (and the later interacting string theory) is of Faraday’s lines 
of force as fundamental objects which become the physical degrees of freedom of the 
electromagnetic field. 

To complement this work on contact interactions in string theory it seems ap¬ 
propriate to return to worldline theories to explore the consequences of allowing 
point particles to interact when their worldlines intersect. Such theories are of sig¬ 
nificant physical interest, since the so-called worldline formalism of quantum field 
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theory [4, 5, 6] expresses physical quantities in a held theory in terms of the quan¬ 
tum mechanics of point particles which trace out worldlines in space-time [7, 8, 9]. 
This technique can also be extended to non-Abelian theories and chiral fermions [10] 
where it may provide insight to the unihcation of the fundamental forces [11]. The 
coupling between matter helds and the gauge held is described in the point particle 
theory as a Wilson-line interaction for the particle worldlines in the presence of a 
background held. For example, the partition function for a single scalar held mini¬ 
mally coupled to an Abelian gauge held, A, is expressed in the worldline formalism 
(which we derive in more detail later in the paper) by an integral over all closed 
curves to 

/ OO jrp p 

— j) [A] ; W [A] = (1.1) 

where S'point [i^] = / dry/cj^ is an action describing the dynamics of a point particle 

and W [.A] is the Wilson loop describing the interaction of the particle with the 
gauge held (we have absorbed the coupling strength into A). The right hand side is 
interpreted as quantum mechanics on the worldline of this particle and it is this hrst 
quantised theory which we propose to modify in this paper. 

Field theory is the conventional framework in which to introduce interactions and 
the local nature of this approach naturally leads to particles interacting upon contact. 
However, the worldline formalism can oher substantial calculational advantages over 
traditional approaches in held theory, especially since it represents a reorganisation 
of the usual perturbative expansion in Feynman diagrams and makes the local gauge 
invariance manifest [12]. It is therefore important to develop worldline techniques 
further and one of the most basic modihcations to the theory must be to introduce 
direct interactions between these particle worldlines. As we shall describe below, 
a modihcation to the worldline theory can be interpreted as inducing a change in 
the underlying held theory, so the results of this program may provide new tools to 
complement the conventional techniques familiar to held theorists. We comment on 
this in section 5. 

Direct inter-particle interactions can be found in many previous publications. 
One of the most well-known approaches is the action at a distance formulation 
of electrodynamics by Feynman and Wheeler [13, 14], originally proposed to ad¬ 
dress the problem of radiation reaction. This built upon earlier work in formulating 
a consistent theory involving direct inter-particle interactions by Tetrode [15] and 
Fokker [16], who described electromagnetic phenomena in terms of interactions be¬ 
tween particles with light-like separation. Ramond generalised this work and found 
a set of consistency constraints limiting the form of the interaction that can be 
introduced into worldline theories [17]. This was further extended to include di¬ 
rect inter-string interactions as well as interactions between particles and strings 
[18, 19, 20, 21, 22]. A number of further theories involving action at a distance have 
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been proposed to describe various other phenomena within this framework - see for 
example [23, 24, 25, 26, 27]. 

The general principle is to couple the worldlines of the particles together by 
adding extra terms to the free particle action, S'point- An illustrative example would 
be to consider a theory of two particles whose worldlines are described by {to) 
and uj^ (rfe) and to introduce [28, 25] 

Sint = 9a9b / dTa / du dja (^a) ' UJb {u) D {uJa - UJb) (1-2) 

Jua JoJb 

as an interaction term in the action^. Here the function D must be symmetric and 
accounts for the relative strength of interaction as a function of particle separation. 
In the literature discussed above D {ua — oJb) has been taken to be supported for 
light-like, time-like and space-like separations, the last of which is the relativistic 
generalisation of an instantaneous interaction. In Feynman-Wheeler theory, for ex¬ 
ample, D {ua — ojb) is taken to be the sum of advanced and retarded Green functions 
of the (space-time) Laplacian. As is now well known, these action at a distance the¬ 
ories are often cast into a form reminiscent of a field theory, although the “fields” are 
not independent variables but are rather defined in terms of the particle trajectories 
and the choice of D. For instance, if we take [29] 

D {uJa - UJb) = -^6 {{Ua - (1.3) 

^ ^ y (Wa - UJb)"^ ^ ' 

which is the time symmetric Green satisfying [d^d^ — m?) D {uJa — uJb) = (uJa — uJb) 
then we may dehne 

(x) = '^9a f dTa UJa i'^a) D {uJa (Tq) - x) . (1.4) 

a 

This satishes the Maxwell equations and gauge condition 

{d^d‘' - m^) (x) = (x); (x) = - ^ f dxa uj^ {xa) 5^ {x - uJa) , 

a 

(x) = 0 (1.5) 


and in terms of A the interaction between the particles takes the form 

^5'a / dTaUJa{ra) ■ A{uJa{Ta)) (1.6) 

a duia 

which shows that the action at a distance formalism contains the same equations of 
motion and interactions as more traditional approaches using field theory. 

^To be precise this provides vector-like interactions between point particles. Scalar interactions 
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can be produced by replacing each by (cof) ^. 
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The proposal we will make will follow the same form as (1.2) except that we shall 
choose D {ua — ujb) = 5^ (cUq — Ub) so as to provide contact interactions between the 
worldlines. This also ensures that, although in principle (1.2) implies the interaction 
is non-local on the worldlines, the particles only communicate when they meet so 
that the theory is local in space-time. In other words we are no longer considering 
action at a distance but we allow particles to interact when they hnd themselves at 
the same space-time position. As stated above, we have previously considered such 
contact interactions between strings, where the theory found application to classical 
electromagnetism and quantum electrodynamics. We now intend to explore the same 
ideas for the case of point particles. 

This article revisits and extends the results of [3] and also generalises that work 
to the case of fermionic particles. It then goes beyond leading order in the coupling 
strength to demonstrate that in fact the full quantum theory of a set of interacting 
point particles is consistent and free of unwanted divergences. We develop the func¬ 
tional approach to one dimensional held theory for consistency with [1, 2] and for the 
generalisation to fermionic particles we will hnd it most natural to form the theory 
in superspace. We will hrst consider a single particle worldline with hxed endpoints 
that is constrained to pass through a given point in space and will relate it to clas¬ 
sical electrostatics and the well-known phenomenon of conhnement. We then repeat 
a similar calculation for spin 1/2 particles to explore the fermionic version of the 
theory before considering an arbitrary set of interacting worldlines. We will see that 
the partition function of this theory is related to propagators of the Klein-Gordon 
operator. 

The structure of this article is as follows. In section 2 the bosonic theory pre¬ 
sented at lowest order in [3] is reviewed before we generalise it to include spin degrees 
of freedom in section 3. In section 4 we also carry out the hrst analysis of the theory 
beyond the classical limit to explore higher order corrections to the result in [3]. Fol¬ 
lowing this a full quantum theory of interacting worldlines is described and quantised 
in section 5. Some supporting calculations on our regularisation scheme are given in 
the Appendix. 


2. Bosonic particles — the classical electric field 

We begin by working in D spatial dimensions and consider a static charged particle 
at position a and an oppositely charged particle at b. The classical electric dipole 
held for this conhguration is the well known solution to Maxwell’s equations in the 
presence of these point particles. In [3] an alternative proposal was made which 
generates the electric held by carrying out an average over a set of curves joining 
the two sources. This concept goes some way to reviving Faraday’s notion of electric 
hux lines and it is this calculation that we now review and extend. 
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An expression satisfying Gauss’ law^ was given in [3]: 



where the integral is taken over any curve C with endpoints at a and b. The form 
of this held is reminiscent of the Dirac string which was introduced to describe the 
held of a magnetic monopole [30, 31] but we shall use it here in the context of 
electrostatics, where C is interpreted as a single line of hux stretched between a and 
b. The expression for E' does not satisfy V x E' = 0, and so cannot represent the 
physical held of a pair of charged particles, but we shall see that its statistical average 
does. This average is over all curves with endpoints hxed at a and b and is dehned 
in reparameterisation invariant form as 

1 /•oj{i)=h / \ 

(D(n;))^ = — / (lj) 6 i / edr — T ] ] S[(jj,e]= / —dr 

Z Ja>(o)=a J Jo 2e 

( 2 . 2 ) 

where D is any reparameterisation invariant functional of the path (we work in Eu¬ 
clidean space). The action is that of Brink, diVecchia and Howe [32] describing the 
dynamics of a bosonic point particle (we take the particle to have zero mass^) whose 
worldline traces out a curve u: (r), which the square of the einbein, e (r), equips with 
a one dimensional metric. The normalisation constant is dehned by (1) = 1. This 
action is invariant under diheomorphisms r —)■ f that preserve the parameter interval 
if under such reparameterisations e^(r) transforms as a metric and n; (r) as a scalar: 

e(f) df = e(r) dr; u}{f) =u>{t). (2.3) 

This idea of averaging a quantity over curves with hxed endpoints is similar in spirit 
to some applications of the worldline approach to the Casimir ehect - in [33, 34] the 
force between two planar bodies induced by quantum huctuations of a scalar held 
was arrived at by counting closed worldlines which intersect both surfaces. Open 
worldlines have also been used to calculate dressed propagators for matter interacting 
with a background gauge held [35]. 

The d-function in (2.2) will appear unfamiliar, both in the context of quantum 
mechanics and the worldline formalism, since it picks out paths of hxed intrinsic 
length^ T. The Coulomb held for the pair of particles was arrived at in [3] by taking 

^We put eo = 1 and denote the electric charge by g so V • E' = qS^ (x — a) — qd^ (x — b) 

mass term could have been included in the action defined in (2.2) via the inclusion of a 
cosmological term fg edr. Usually it is necessary to do so in order to remove a divergence in 
the functional determinant from the integral over a; by a renormalisation of this bare mass ttiq to a 
physical mass m. In the current case the divergence is cancelled by 2 and since E' does not involve 
the einbein the effect of including m would also be cancelled by the normalisation. 

^In the worldline approach it would be more common to interpret T as a modulus which remains 
after gauge fixing the einbein. As has been discussed in [3], the average defined in (2.2) retains 
reparameterisation invariance, despite the d-function restricting the functional integral to the subset 
of einbeins which give the paths a fixed length. 
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the average of E' in the limit as T —)■ cx): 


^im (E'(x))^ = ^im ^ J edr — T^ j dr h; 5^ (n; (r) — x) e 

(2.4) 

Taking this limit means that the lines of flux can fluctuate on a macroscopic scale. 
This is necessary because the held at the point x is constructed by adding together 
the contributions from all curves that start and end at the location of the charges 
and which intersect the spatial point in question. In order that the held permeate all 
of space, rather than being conhned to regions that are of a comparable scale to the 
distance between the points (in comparison to \/T), we must allow the worldlines 
unconstrained huctuations. For comparison to [1, 2], which involved averages based 
on Polyakov’s formulation of spinning strings, T plays a similar role to the string 
tension, a'. In our earlier work we found it necessary to take the limit in which 
the strings were tensionless {a' —>■ oo) which allowed the huctuations of the string 
worldsheet to be macroscopically large. This ensured that the averages were not 
dominated by the minimal spanning surface of the string boundaries and allowed an 
Abelian gauge held to be supported throughout space-time. In the current context, 
the large T limit plays an analogous role, ensuring that the held lines have access to 
all of space. In section 4 we will return to discuss the consequences of hnite values 
of T, exploring also the low T limit for both bosonic and fermionic particles. 

In the current context of point particle worldlines, it is also possible to interpret 
T as an ehective temperature so that the functional average dehned by (2.2) can be 
regarded as a thermal statistical average. This was explained in [3] and can be seen by 
exploiting the reparameterisation invariance to set the einbein to be constant along 
the worldline. The constraint then imposes e = T and the action can be written as 

where ^ = T~^ and l-L = uj^/2dT. Then integrating over all curves with hxed 
endpoints can be seen as allowing thermal huctuations of the worldlines, the relative 
weighting of each conhguration being related to the temperature. As discussed above, 
desiring that large-scale huctuations are not penalised in the average requires the 
limit T ^ oo. We will return to this interpretation at the end of this section when 
we consider the ehect of applying our theory in Minkowski space via a Wick rotation 
of one of the spatial coordinates. In the meantime we think of the constraint in the 
geometrical context of holding the length measured along the worldline constant and 
now return to demonstrate that the average of E' results in the Coulomb held. 

In [3] this result was shown using techniques derived from canonical quantisation. 
Here we conhrm the result using the functional methods we will use for the fermionic 
generalisation of this claim. To do so we must address the overcounting caused by 
the reparameterisation symmetry. Following Polaykov [36] (see also [37, 2]) we recall 
that for open curves the measure on the space of one dimensional metrics can be 
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written _ 

dr e 1 Det {DW ); DV = ^ {Ve) (2.5) 

J dr 

where c represents a scaling of the einbein and V generates an inhnitesimal repa- 
rameterisation. This volume element follows because any metric can be written as 
a combined scaling, generated by c, plus reparameterisation, generated by V, about 
some reference metric, e. We choose to expand about e = T, a constant, whereby 
DWV = —(^^)) = The inhnite product over the eigenval¬ 

ues of this operator can be (^-function regulated as in Appendix A, and the volume 
element becomes 

= dT^V, ( 2 . 6 ) 

whilst the constraint on path lengths becomes 5 (T — T). Since the components of the 
functional average are taken to be invariant under reparameterisations the functional 
integral with respect to V evaluates to the volume of the reparameterisation group, 
which cancels with the corresponding contribution from the normalisation constant. 
So for the average of E' (x) we must determine 


= dc^V 



{I(x)) = |, j 


d^k 



dri 


(2.7) 


where Z' is what remains after the volume of the reparameterisation group has been 
cancelled from Z and we have used the Fourier decomposition of the ^-function. The 
insertion that arises has a familiar form - it is the one dimensional version of the 
vertex operator used in bosonic string theory and it frequently appears in calculations 
in the worldline formalism: 




t) = 


T] e 


ik-io{T) 


(2.8) 


Unlike in string theory we integrate this operator over all momenta. From the point 
of view of the one dimensional quantum theory (2.7) is the amplitude for the path 
from a to b to pass through the point x. Before calculating this expectation value we 
also note that the structure of the insertion allows us to constrain its dependence on 
momentum - if we contract (2.8) with and integrate over r we hnd a contribution 
only from the endpoints of the domain 





i - e- 


(2.9) 


which we shall refer to as the generalised Gauss’ law. 

The insertion can be generated by introducing a source, j (r), and dehning 
J (r) = —q^ — ikd (r — ti). Then the above equation becomes 


1 



d^k 


dT^ e 


— 2kx 


(^l) 


3 /o 2T +^' 


j=o 


( 2 . 10 ) 
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We split ui into its classical part in the absence of a source and a piece which absorbs 
the source and accounts for the quantum fluctuations n; (r) = cjc (r) + d; (r). Here 
LO satisfies the source-free classical equation of motion = 0 with endpoints at 

a and b: 

cj (r) = a -|- (b — a) r. (2-11) 

The Dirichlet boundary conditions mean that Gj (r) is required to vanish at the 
endpoints. It can be split into a classical piece, cdc, satisfying = J (t) and 

a quantum fluctuation, lj (r). Integrating over d; leads to a functional determinant 
(which we evaluate with (^-function regularisation in Appendix A) and because the 
path is open there are also boundary contributions. We find 


(27rr) 


_D (b-a)^ 
2 0 2T 




d^k 


dn ^ 


>0 


(h) 


n 


d 


exp --k^G{Ti,Ti) -iq drk-j (r) 


'0 


dr 


r/r 


<ir<;r'j(r).j(r')|;TG(r.r')) 


exp ( -q / drj (r) ■ (r) + ik ■ (n) 


. ( 2 . 12 ) 


j=o 


In the above equation G(ti,T 2 ) is the Green function on the interval [0,1] with 
Dirichlet boundary conditions G (0, r 2 ) = 0 = G (1, r 2 ). Its explicit form is 


T 

G (ri, r2) = -— {\ti - T2\ - {ti + T 2 ) + 2tiT2) ; 


1 d'^G 


d{ri-T2). (2.13) 


and G{ti,ti) is the coincident limit, which in one dimension is finite. Worldline 
Green functions play an important role in first quantised formulations of field theory 
and as such have been analysed extensively [38]. The numerator of the prefactor 
in (2.12) is the heat-kernel for the free particle and this cancels exactly with Z'. 
Garrying out the functional differentiation and setting j = 0 yields for the average 
in momentum representation 


(I(k)) = -g 




l,kAG(r..r.) 


g-|DG(Tl,Tl)gik-tJc(Tl) 


(2.14) 


Since we will eventually take the limit T —)■ cx) it is useful at this point to extract the 
T dependence in order to make an expansion in powers of Dehne then G (ri, T 2 ) = 
^G (ti, T 2 ) and note that 


Giji.Ti) = -Ti (ri - 1). 


(2.15) 











This function vanishes on the boundary and is increasing (decreasing) for r < 1/2 
(r > 1/2). In this parameterisation (2.14) becomes 



. /X T d ~ , . 

/ dri 

(t) - —G (ti, Ti) 

Jo 

2 dTi 


g-iA;2TG(ri,Ti)gik-a)c(Ti) 


( 2 . 16 ) 


In the large T limit the suppression caused by the exponent exp y—^k^TG (ri,ri) 
causes the contributions to the integral to arise primarily when ri —)■ 0 and ri —)■ 1 
(where the coincident Green function vanishes). We therefore expand the field lo 
about these points and integrate a small distance, h, along the worldline. At lowest 
order in we have: 


- q 


- q 


dvi 


dri 


'l-h 


^c(O) - ^fk^G(ri,ri) 
2 CLTi 

^c(l) - '^ik^G{Ti,Ti) 
2 CLTi 


3-1 A:2rG(Ti ,Ti) gik-tJc (0) 


(2.17) 


and the damping caused by the form of the coincident Green function in the exponent 
allows the integration regions to be extended by setting h = | as we take T —)■ 
00 . Each integral has two terms, the second of which provides the leading order 
contribution: 


Tik 


CLTi 




1 


ik 


( 2 . 18 ) 


where the last line holds as k^T —)■ cx). Noting that for 0 < r < 1/2 we have 
G {t,t) > I the first terms of each line in (2.17) can be bounded 


rG(ri,ri) < / 


-Ifc^Tri 


k^T L 


-ik^T 


e 8 


(2.19) 


which is O Putting this together with the contribution from the other end of 

the path we arrive at the momentum space expression 


(I(k)) = ^(e*^"-e‘*^-^)+(p(^ 


FT 


( 2 . 20 ) 


We check our answer by compatibility with the generalised Gauss’ law: dotting with 
k produces the expected contribution (2.9). Finally, taking the limit T —)■ 00 and 
setting D = 3 the result can be written in position space as 


(I(x)> = Xv 


X — a 


+ 


lx — b| 


( 2 . 21 ) 


9 

















which is indeed the classical dipole electric held. This average therefore determines 
Foi for static oppositely charged particles. In [3] the case of magnetostatics was also 
considered for a hxed closed current carrying wire, and the time varying situation 
was also examined. These cases require treating the curve C as dynamical so that 
the natural weight becomes not the action of a point particle but that of extended 
objects, naturally leading to the use of string theory. This picture was explored 
further in [1, 2] where a theory of tensionless spinning strings interacting on contact 
was shown to be equivalent to the Abelian theory of quantum electrodynamics. 

The calculation presented above provides an interesting method of determining 
the classical static dipole electric held, albeit somewhat unconventional. The utility 
of the functional approach is the ease with which it can be extended. Before gen¬ 
eralising to fermionic particles we show that a simple change can instead generate 
the classical electric held due to a static point particle. This is desirable since it 
is presumably more useful to deal with a single particle rather than be constrained 
to dealing with oppositely charged pairs (except when considering the worldlines of 
virtual particle / anti-particle pairs). For a single particle at the point a we proceed 
as above with the exception that we constrain only one end of the worldlines. This 
is a simple change of the boundary condition at the upper end of the interval; the 
variation of the action shows that the only other consistent choice we can make is 
the Neumann condition ch (1) = 0. 

There are three effects of this change. Firstly the determinant of the kinetic 
operator becomes independent of T (see Appendix A). This has no relevance 

because it is cancelled by the same change in Z'. Of more importance is the change 
in the Green function. With the new boundary conditions we hnd 

G' (a, 7 - 2 ) = (|ri - r2| - (ti + T 2 )) ; G' (ri, ri) = Tti. (2.22) 

Note that the coincident Green function now only vanishes at ri = 0, the end of 
the curve tied to the particle at a. Finally the classical solution must clearly differ; 
now LJc = a satishes the equation of motion and boundary conditions. In particular 
cj = 0 and there are no boundary contributions from the classical action so that 
(2.16) becomes 

(I(k)) = q (2.23) 

Jo ^ 

Dehning G'(ri,r 2 ) = TG'(ri,ri) we could again consider the above equation for 
large T whereby the integrand is suppressed by the exponent exp |fc^TG' (ri, ri) j 
except for ti ~ 0. The integrand is a total derivative however so we can determine 
the exact answer. We obtain 

(I(k)) = |^(l-e-5'=^^)e**^'" (2.24) 
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which in the limit T —)■ cx) has inverse Fourier transform equal to the static electric 
held of a point particle of charge q at the spatial point a. 


= (2.25) 

dvr |x - a| 

It is interesting to note that (2.24) gives the exact average at hnite T. We will come 
to explore corrections to the classical helds for both of the above conhgurations of 
particles below. We have already commented on the calculation of the static magnetic 
held of a closed loop of current carrying wire given in [3]. This required a functional 
average over surfaces with boundary hxed to the wire and already invoked the use of 
string theory. There the two dimensional worldsheet Green function was required to 
vanish at both ends of the string since the endpoints were constrained to lie on the 
boundary. The analogous method of applying Dirichlet boundary conditions to only 
one end corresponds to opening up the wire into a small segment and ensures that 
the only contribution to the average comes from a strip close to the end of the string 
hxed to the wire. The average then yields the Biot-Savart law for that segment of 
wire. Mixed boundary conditions are discussed further in the full interacting string 
theory in Appendix A of [2] where it is shown to provide propagators of the held 
theory. 

For the most part we took D to be arbitrary, only specialising to D = 3 spatial 
dimensions for the sake of compatibility with [3] and an illustration of some of the 
physical content of the average. In a four dimensional space-time we will have to 
deal with 4 bosonic coordinates so we append to the three helds u considered 
above. The Euclidean average is then constructed over all paths which intersect the 
spatial point = (a;°,x) whose endpoints are hxed to the particles at = (a°,a) 
and = (6°, b) . The average takes the form 

(x)) = (q j (^ (^) ■ ( 2 - 26 ) 

If we restrict attention to the case of a single point particle at position discussed 
above the generalisation of (2.24) provides the Euclidean space average 

{I^^ (k)) = ^ (^1 - (2.27) 

We can interpret this in Minkowski space by treating the boundary data as some 
hxed point on the worldline of the particle and Wick rotating the above result. In 
the T —)■ CX) limit then 

= -i2— (2.28) 

\x - a\ +te 

where we use the ie procedure to specify the positions of the poles. 

The physical interpretation of this result is less obvious because it relies on 
the choice of (also had we included a second particle). If we return to a static 
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picture then consider an observer at in the rest frame of the charged particle. The 
calculation in D = 3 spatial dimensions was an average over all possible particle paths 
starting at a and passing through x but in four dimensional space-time they are also 
required to pass though x at the time = t, say. We are restricting our attention to a 
static conhguration and suppose that we ought to integrate over all possible starting 
points on the worldline at which the path could begin whilst still passing 
through x^. Recall that the parameter T can be interpreted as a temperature, which 
in Euclidean space implies that our calculation is a thermal average. In section 4 
of [3] an argument was given based on the construction of thermal Green functions 
that the retarded solutions to Maxwell’s equations are inherited from the Feynman 
propagator if the Minkowski space calculation is seen as a hnite temperature quantum 
expectation value. Application of this procedure to the current problem provides 


(^)) 


q ^ 5 (x° — a° — |x — a|) 

~A ^ ^ 

Att X — a 


(2.29) 


so that the only contribution comes from paths whose endpoints are joined along the 
causal light-cone. In the rest frame of the particle its worldline has constant a and 
integrating (2.29) with respect to a° yields the static electric held expected at x^ = t: 


Fqo — {Iq (x)) — 0; Foi — {li {x)) — ^ (2.30) 

dTT |x — a| 

The same resnlt could be arrived at by integrating (2.28) over a° immediately if 
the contour for the a° integral is chosen to fall above the poles on the real axis; 
the discussion of the thermal average in [3] can thus be seen as justihcation for this 
choice of contour. The result (2.30) is then familiar as the held created at position 
(x°,x) due to a point electric charge which had position a a time |x — a| earlier, 
which appears in the more familiar action at a distance theories discussed in the 
introduction. 

It should be stressed that this procedure yields the correct static held but is 
not applicable for the general time dependent problem. This has been dealt with in 
[3] and [1]. It requires string theory to correctly describe the dynamics of extended 
curves whose boundaries are hxed to the worldlines of the charged particles. The 
contact interaction is then between fundamental strings, between worldsheets rather 
than worldlines. The static problem considered above is an unusnal way to arrive 
at the electric held but is nonetheless of interest because of its straightforward gen¬ 
eralisation. In the following sections we hrst include spin degrees of freedom before 
returning to an analysis of the result when the parameter T is taken to be hnite. 
We then generalise the work to form a full quantum theory of point particles with 
contact interactions. 

Following on from the bosonic theory we wish to provide an extension to the 
theory in [3]. The following sections comprise the new contribntion in this article. 
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We have three aims in sight. The hrst is to extend the work to fermionic helds and 
the second is to determine the corrections to both results which are present at finite 
T. We hnally ask whether the interaction exponentiates as in [1,2]- we shall address 
this issue in section 5. 

3. Fermionic particles 

In this section we consider a theory of spin 1/2 particles and generalise the theory 
presented above for application to this case. Spinning particles have been considered 
in action at a distance theories in the past [39, 40, 41] and it is natural to expect 
some spin-dependent modihcations to the results found in the previous section. We 
continue to work with massless particles for simplicity and it will also prove most 
convenient to work in four dimensional Euclidean space (this is more natural for 
fermionic theories - we discuss how we may relate this to the three dimensional case 
below). To deal with fermions it is useful to construct the theory in superspace. 
We therefore introduce the Grassmann variable 6 to extend the parameter domain 
r —)■ (r, 6). We further introduce the scalar superfield® 

X{r,9) = u {t) + 9e^ (3.1) 

and the supereinbein 

E (r, 0) = e (r) + (r) X (r), (3.2) 

where ijj is the superpartner to u and y is the gravitino. We also define the su¬ 
perderivative 

D = de + 9dr. (3.3) 

Under the local supersymmetry transformations parameterised by V (r), the gener¬ 
ator of reparameterisations, and rj (r), a Grassmann function generating pure super- 
symmetry transformations, 

T ^ T + V (t) + 9r/ (t) ; 9 9 + 7] (r) + ^9V (r) (3.4) 

the superderivative transforms homogeneously 

DX ^ A (r, 0) DX (3.5) 

and the supereinbein transforms as 

E^A2(r,0)E (3.6) 

®The strange looking factors of es are necessary for the action given in the text to reduce to 
the familiar action of Brink, Di-Veccia and Howe [32]. Another convention for the superhelds exists 
where we scale ip ^ and x —>■ 
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where A (r, 6^) = 1 + (r) + 9f] (r). Requiring the integration measure to transform 

as drclO —)■ (r, 9) dTd9 the following action is invariant: 


5[E,X] = - 


dTd9 


(3.7) 


Integrating over 9 allows this to be cast in the more familiar component form® 

^ j dr + tp ■'tp (3.8) 
z j e 

Under canonical quantisation the equations of motion for the auxiliary helds y and 
e lead to the hrst class constraints p ■ p) = Q and then = 0 respectively. Here 
p^^ = Oj^ je is the momentum corresponding to the held uj. The fundamental anti¬ 
commutation relations ip'^} = 5^'^ can be solved by taking ^p^ = ■^'y^ so that on 
the state space the former of the constraints enforces the Dirac equation on physical 
states and the latter informs us the particle has zero mass [32] . Below we shall pursue 
again the functional quantisation of this action. 

We also require the supersymmetric generalisation of the interaction and the 
constraint on the intrinsic length of the path. The natural invariant interaction term 


IS 


I(i) =« / dTdeDXS^(X-x), 


(3.9) 


whose functional average will be taken over all worldlines with endpoints attached 
to hxed charges. Fourier decomposing the 5—function and integrating over 9 puts 
this into the form 

9 / / <iT (lO - ei,ik ■ V.) (3.10) 

which is analogous to the supersymmetric vertex operator familiar in the context of 
the spinning string and the worldline formulation of spinor held theory: 


l/' (t) = (r) - e (t) i/.i' (t) ik ■ i/> (t)) 


(3.11) 


Such an expression was examined in a string setting in [1, 2] but for now we continue 
to explore the point particle theory. We contend that the generalisation of the electric 
held part of the held strength tensor should be generated by a functional average of 
(3.10) with an appropriate weight. This weight will be of conrse that corresponding 
to the action in (3.8). Note also that the anti-commuting nature of ip ensures that 
the generalised Gauss’ law (2.9) still holds for this interaction. 

The hnal element we need is the supersymmetric version of the constraint on path 
lengths. Previously we inserted edr —T) into the functional average bnt this 

®The supersymmetry transformations of the component fields follow from those of the superhelds 
and (in the absence of reparameterisations) are: ^rjip = f ~ = VW ^vX = 

2f]. Under reparameterisations, p \s a worldline scalar like w, whilst y transforms as e. 
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changes under supersymmetry transformations. Introducing an arbitrary Grassmann 
number S, the natural invariant quantity is 

6 (^j dTd9E'2 (3.12) 

which imposes a Grassmannian constraint on y rather than on the metric. In the 
massless case considered here it is not possible to construct a local function of e and 
y which is supersymmetric. Instead we follow Polyakov [36] and give a superspace 
version of the non-local and invariant quantity he termed the superlength; 

dTidOi (^Ti, 9i) D 1 D 2 G (ti, 9i] T 2 , ^ 2 ) (r2, ^ 2 ) dT2d92 — , (3.13) 



where G (ri, 6 ^ 1 ; r 2 , 6 * 2 ) = |ti — r 2 — 6 'i 6 * 2 | is a superinvariant Green function and Di 
is the super-derivative acting on the parameters {Ti,9i). In components the first of 
these constraints takes the form 




drx{r) - 2 ^ 


whilst the latter can be written 


'^TeiT)- - / dTiX (ti) sg (t, - T2) X (ti) dr^ - T 


(3,14) 


(3,15) 


with sg (r) = equal to the sign of its argument. We shall require the superlength 
constraint in the functional average because it will be seen to provide the appropriate 
hxing of the einbein and we will also impose the complementary constraint (3.14) 
which will similarly fix y. 

To continue the calculation it is necessary to determine the measure on the space 
of the gravitino and ijj and also to specify the boundary conditions we will use. It is 
well known that the purpose of the helds is to represent the 7 matrices (in canonical 
quantisation we have seen it suffices to take oc 7 ^), which is why it is desirable to 
work in a four dimensional space-time. For instance, consider the partition function 
of a free fermion -0 coupled to a source C,. Integrating the exponential of the action 
over 0 yields an object with spinor indices and choosing the boundary conditions 
on the integral allows us to extract each component of the result. Specifically in 
Appendix G of [ 2 ] we have shown that, for example. 



'i/>2=—*7 i;V’4=*73 


72=171 ;'i/> 4 =—*73 
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(3.16) 


where ^ represents time ordering along the worldline. We also showed that the 
volume element for y can be written 


= dxo^rj 


_ 1 

Det “2 



(3.17) 
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where Xo is the constant piece of x reqnired for consistency with the bonndary con¬ 
ditions of open worldlines. This volnme elements follows becanse y can be expressed 
as a snpersymmetry transformation, generated by rj, plus a change proportional to 
e about some reference, y. We gauge £x by expanding about e = T - a constant ~ 
and y = yo- 

These are consistent choices which cover the physically distinct conhgurations 
and on this gauge slice the action becomes 

^ J + - ■ -0 (3.18) 

whilst the interaction vertex is given by 

(cj^ (r) — (r) ik ■ ip (r)) (3.19) 

Similarly, using (^-function regularisation for the functional determinants (see Ap¬ 
pendix A), the volume elements become 

= dTdxo^V^V- (3.20) 

The super-length constraint reduces to 

d(T-T), (3.21) 

and the analogous version for y becomes 

d(yo-S), (3.22) 


which can be used to carry out the integrals over T and yo. We shall again eventually 
take T to infinity so we can expand in powers of ^ and we will take the dimensionful 
Grassmann parameter S to vanish. The action and insertions are locally supersym¬ 
metric so that the integrals with respect to V and rj will evaluate to the volumes 
of the reparameterisation group and supersymmetry group respectively. These con¬ 
stants, albeit formally divergent, are cancelled by their counterparts if we normalise 
against the bare partition function. 

We pause here to derive an important result using (3.16). We consider the 
expectation value over the fermionic helds {^jJ^ (r) "0^" (r)) where we have attached 
the spinor indices a and 0 to the beginning and end of the worldline respectively^. 
The insertions can be generated by introducing a fermionic source rj and carrying 
out functional differentiation, after which we set rj = to produce the linear 

term in the action. Then (3.16) gives 


{V (r) 0^ (r))„^ ^ 


5 


5 


6ri>^ (t) 67]^ (t) 


^ ('g/rfr^.7 


a /3 


' 2\/2T 


(3.23) 


^We have written the fields at equal times but there is of course a time-ordering issue. We take 
the convention that the field to the left is understood to be evaluated at an infinitesimally greater 
time, e, than that on the right, after which we take the limit e —>■ 0"^. 
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where the constant Z is determined in the appendix. If we impose S = 0 the resnlt 
of the fnnctional differentiation and the integral over xo is 

W-'‘(r)r(T)) = ^(5'“'-7'‘7") (3.24) 

which crncially is independent of the held u, since by setting S = 0 we have decon- 
pled the helds u and "0 in the action. This resnlt is consistent with the symmetry 
properties of the Grassman helds and will play a key role in the forthcoming calcn- 
lations. 

We now retnrn to (3.10) and take its fnnctional average over worldlines hxed to 
charges at and b^. Carrying ont the integral over T we are left with the gange 
hxed Fonrier space expectation valne 

(F {k)) = ^ I (xo ■ 0 ) 

dn - y - Yk ■ 7 )^ -fo (3.25) 

where for the last line we have integrated over ip to prodnce the gamma matrices and 
integrated over xo before pntting S = 0. It remains to determine the expectation 
valnes of these integrands and carry ont the integral over ri. The integral over cj can 
be fonnd to give 



q [ dn (ri) -—{k^^- -f>^k ■ 7 ) - Uk^^G (n, Ti)] 

Jo V 2 2 OTi J 

(3.26) 

Making fnrther nse of the redehnition G (ti, T 2 ) = TG (ti, T 2 ) this can be cast in the 
form 


Ti 


d 


q drA (ri)- {k^^ - -f^^k ■ 7 ) - -ik>^—G {n, n) TG{T^,r^)^ik-Mri) 

Jo V 2 2 OTi J 

(3.27) 

From the exponent exp (ti, ti) j it is easy to nnderstand that the latter 

two terms in brackets will contribnte at leading order in ^ and that the hrst term 
will once again be sub-leading. In the limit of large T we can follow the calculations 
of the bosonic theory and approximate the integral by evaluating the classical path 
Uc on each boundary and then integrating a short distance, h, along the worldline. 
We calculate 


'•h 


'0 


Tt 


T 


d 


dri ( ( 0 ) Yk ■ 7 ) " 2 


1 

A-h 


Ti 

y 


+q / driicuY^) — tt (k^ ~ Yk ■ x) ~ Yk^^ —G(ri,ri) e 2 


T 

Y 




^-\k‘^TG(Ti,Ti) ik-0Jc{l) 


(3.28) 
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and as before the suppression caused by the exponent at large T allows us to extend 
the integrands by setting h = 1/2. Carrying out the integrals in the large T limit 
leads to 


{I^{k)) = q 


i 


'-y^k ■ 7 



- + O 



(3.29) 


In the limit T —)• oo only the first term contributes and so we have determined 

= (3.30) 

We may double check our work by contracting with to verify against (2.9) that 

the index structure is correct. In position space the expression above becomes 

= (3.31) 

which we take as the generalisation of the Coulomb held found in the previous section. 

We conclude this section by commenting on how this approach could be applied 
to a three dimensional space. The main obstacle to this case is the structure of 
(3.16) which imposes boundary conditions relating pairs of the Furthermore 

the proof of this equation presented in [2] does not hold for the three dimensional 
version of the Clifford algebra because the chirality operator does not anti-commute 
with the other matrices; in particular the trace of a product of an odd number of 
gamma matrices does not in general vanish. We present two alternatives to address 
this issue. The hrst is to use the symmetry of the problem to choose a coordinate 
system such that the two charges are placed on the z = 0 plane and to restrict our 
attention to determine the field only on this plane. Then we need consider a super- 
field, Lo, which consists of only two components and the fields ip essentially become 
the a-matrices cxi and a 2 '- 

(r (k)> = ^ I ^ dn (cj* - ■ ^P) e*k-e-/o 

^ dn - y (fy - fyk ■ cr)^ 

(3.32) 



where i G {1,2}. Carrying out the integral over the two-dimensional field n; and 
following the same steps as the four dimensional case leads to 

(I (k)) = ^cr k ■ cr . (3.33) 

fv 

It is more satisfactory to instead modify the four dimensional theory in such a 
way that statements can be made for the three dimensional case. This can be done 
in two equivalent ways. The three dimensional version of the vertex operator is 

Vfc = (ch — e^pik ■ 'ip) (3.34) 
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We could use this as an insertion in the four dimensional theory by dehning in 
some inertial frame = (0, 1) where the 3-vector I is constructed out of the three 
dimensional vertex operator above. Then (3.16) continues to hold except that the 
part of the expectation value which does not cancel with Z contains only a three 
dimensional source so the right hand side of that equation becomes 

^ (3.35) 


This is equivalent to beginning with an entirely three dimensional theory but intro¬ 
ducing a further pair of helds Uq, Into the integral of some functional of the three 
dimensional fields we introduce a supersymmetric invariant factor in the numerator 
and denominator as follows: 


1 


j (t^, ' 0 ) e 


fo dr ^ + 


f ^CjQ^'ipQe /o 2 T + 2 '’^ 0 'ho 

(3.36) 


We then combine the denominator of the fraction with the three dimensional nor¬ 
malisation to form 


J UjS^'lpd'XQ e /o 2T + 2 2t‘^'’’^ 2t‘^0-V’0 


(3.37) 


where the integration is now over four bosonic and four fermionic helds. We apply 
the boundary conditions for open paths to both the numerator and denominator of 
(3.36). The integration in (3.37) produces an expression dependent on xo which feeds 
back into (3.36). Carrying out the integral over xo the four integrations with 
respect to "0 produces two terms which conspire to yield 

^ j n (u;, 7 ) e“Io(3.38) 


where the normalisation Z' is precisely the correct factor to ensure that (1) = 1. 
This is what we would calculate if we were to take the expectation of = (0,1) 
as dehned above in a theory with four pairs of helds. It is easy to modify (3.30) to 
respect this change: 

(I (k)) = II 7 k ■ 7 (3.39) 

K 


which is to be integrated with respect to the three-vector k. 
representation 


7 ‘ 



We may choose the 


(3.40) 


to show that this approach involves the Pauli matrices in a similar way to (3.33), but 
now we need four-index spinors. With this representation 7 k ■ 7 is block diagonal 
with the two blocks taking the same form as (3.33), dihering from one another by a 
sign. 
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4. Analysis at finite T 


The classical fields of the previous sections were found in the limit that the dimen¬ 
sionful parameter T was taken to be large compared to momenta. We discussed 
above how this limit can be interpreted as making the intrinsic length of the world¬ 


lines stretched between the charged particles macroscopically large. An interesting 


question is to ask about the form of the statistical average for finite T in order to 
explore how the fields change. In this section we give the subleading correction to 
the fields at large T and also consider the opposing limit T —)■ 0. We shall do so 
for the lowest order interaction which generates the classical fields. The geometri¬ 
cal interpretation of changing T is to explore the contribution to the average from 
worldlines of different fixed lengths, but as we have discussed earlier it may also be 
considered as varying the temperature at which the thermal average is calculated. 
Therefore, at lower values of T, we expect that large fluctuations of the worldlines 
will be penalised and so the functional averages will become dominated by curves 
close to the shortest path between the charges. 

4.1 Corrections in the bosonic case 

We first analyse the bosonic particle with mixed boundary conditions whose high T 
limit produced the classical electric field of a point particle (2.24). That equation was 
exact in T and its position space form is found by carrying out the inverse Fourier 
transform. We shall specialise here to the three physical spatial dimensions. Aligning 
the z-axis parallel to the spatial separation x — a the angular integrals of (2.24) can 
be carried out to produce 



(4.1) 


At large T this can be expanded in powers of ^ and we find 



X — a 


1 



(4.2) 


valid for |x — a|^/T <C 1. The leading order correction arises from the third term in 
brackets: 



giving the finite T deviation from the inverse square law. 


20 














The other limit of interest is at small T so we consider an expansion of (4.1) 
about T = 0. The change in functional form is more dramatic because we hnd 


/!( \\ - -1x7 1 / ^ 

(I (x))j, — V 


1 2 


2 I lx 


X — a 


+ o{t 


(4.4) 


which carries a different power of the spatial separation. Carrying out the differen¬ 
tiation we acquire 


47r \ X — a 


4T 2 (x — a) / 5 

^^ + 0 (T5 

X — a 


(4.5) 


in the limit that |x — a|^ /T S> 1. 

The corrections to the dipole held are less trivial to determine because of the 
time dependence in uj,,- For large T the integrand is still dominated by contributions 
at either end of the curve. We hrst expand about r = 0: 

e^k-^^dn) _ g*k.a /i + ■ (b - a) n - y (k ■ (b - a))^ rf + .. .V (4.6) 


When integrated against exp [^—^k'^TG{Ti^Ti)) each extra power of ti results in a 
further power of The hrst term in the square brackets of (2.16) is already 

subleading in 4 so its contribution to the leading order correction comes from the 
hrst term in (4.6). Integrating this from the boundary to Ti = | gives 

plus corrections exponentially suppressed at large k^T. At hrst order the second term 
of (2.16) provided the classical dipole held and its leading order correction arises from 
the O (ti) term in (4.6). Straightforward integration of this evaluates to 

+ o (r-^) (4,8) 

up to further terms which are again exponentially suppressed. It is instructive to 
combine these two terms as follows: the ith component is given by 


2q rw 
Y F 


lY 


{b - ay e^k *^ 


(4.9) 


which highlights the transverse nature of the correction. The index structure of (4.8) 
means the integral with respect to k has the form 


r Yk 4k.(a-x) 

J (27r)^ k^ 


+ B 


(a — x)* (a — xY 


(4.10) 
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for some constants A and B of dimension [length]”^. They can be determined by con¬ 
tracting each side of the above eqnation hrst with dij and also with (a — x) • (a — x)j. 
The hrst of these gives 


3A + B 


d^k _ -1 

( 27 r)^ dvr |x — a 


(4.11) 


and the second yields 


a — xf {A -I- B) 


(fk (k ■ (a - x)f 

(2„f 


(4.12) 


Choosing the z-axis to align with x — a again allows the angular integrals to be 
done and we hnd that the remaining integral with respect to the magnitude of k is 
proportional to 


1 [°° sin(/c|x —a|) 2 cos (/c |x — a|) 2sin(fc|x —a|) 

(27r)^ Jo . |x — a| (fc|x —a|)^ (fc|x —a|)^ 


(4.13) 


Integrating the second term in square brackets by parts once and the hnal term by 
parts twice serves to cancel the hrst term. It is then easy to check that the integral 
evaluates to zero so A -|- i? = 0 and (4.10) evaluates to 


-1 

Svr lx — a 


Sij _ 



(4.14) 


So (4.8) evaluates to 


Q 


AttT |x — a| 


(b-a) 


X - a 


(x — a) ■ (b — a) 


X — a 


o ( t - 2 ) , 


(4.15) 


to which it remains to add the contribution with the analogous calculation for the 
other end of the curve hxed to the point b. The hnal correction at order ^ is 
determined to be 


Q 

AttT |x — a| 
Q 

AttT |x — b 


(b — a) + (x — a) 
(b - a) + (x - b) 


(x — a) ■ (b — a) 

I |2 

|x — a| 

(x — b) ■ (b — a) 
I iTTs 


(4.16) 


and is easy to check that this is divergence free. This is the variation from the classical 
dipole held for large values of T which is present in this model. It is interesting to 
note that with the relaxation of the limit to large but hnite T comes dependence on 
the direction (b — a) which is independent of the spatial point in question. In Fig 1 
we provide an example of the deviation from the well-known dipole held for hnite T 
by plotting the streamlines of the electric held. 
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Figure 1: The field lines for large but finite T demonstrating the small deviation 

2 

from the classical dipole field. In this plot we set T = 10 |b — a| . 


We finally turn to the low T limit of the dipole field for which it is more convenient 
to carry out the integral with respect to k of (2.16), before looking at an expansion 
in powers of T. Here we shall see more striking dependence on the direction of 
separation between the two charges, since as T —)■ 0 one expects the fluctuations of 
the worldlines to be almost completely frozen out. Then the field will be concentrated 
in a small region (measured with respect to T) about the straight line from a to b. 
The result of carrying out the inverse Fourier transform of (2.16) is 


(I ix))T = 


dri 


ClnyJo (^TG(ri,ri) 


G (ri,ri) 
2G (ri,ri) 


X - UJr.) - u;. 


X - LJr 


exp 


2rG(Ti,Ti) r 


(4,17) 

Examining the T-dependence of this expression we see that the limit T —)■ 0 provides 
a representation of the 6 function so that the field is supported only on the classical 
straight line path joining a to b; 

6^ ( [ cLti (x - Uc (ti)) 


= 5" 


dTi (x - (a + (b - a) n)) 


(4.18) 


This makes sense since the T —)■ 0 limit implies the worldline’s intrinsic length is 
small, so it must find the configuration that minimises the intrinsic distance traced 
out by the particle. Alternatively, as we have stated above, there can be no thermal 
fluctuation about the classical solution to the equations of motion. 

To determine the form of the field at finite T we shall employ a Laplace approxi¬ 
mation. For small T the contribution to the integrand of (4.17) is concentrated about 
the positions of the maxima of the exponent exp ^ precise version 
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of Laplace’s method we shall invoke is that for small T and arbitrary well-behaved 
functions / (r) and g (r) 


^ drf (r) = (tq) (1 + 0 (T)) (4.19) 


where the Tq G [0,1] are determined by the condition that g (tq) be a maximum. In 
this case a straightforward calculation shows that the exponent of (4.17) attains a 
single maximum within the integration range at® 


To = 


|x — a| 

X — a| + lx — b 


(4.20) 


If the spatial point x lies on the line joining a to b then the exponent vanishes 
at this value of Ti, in agreement with the T —> 0 limit which provides ^-function 
support on this line. Furthermore the first term in square brackets of (4.17) vanishes 
because x — lJc ("To) = 0. In this case the held is approximated at lowest order in T 
by 


(b - a) 


27rTG (ro, tq) |b - 


<1 

b- 

- a 


2ttT 

X - 


X ■ 

-b| 


(b — a) 


(4.21) 


which has its minimum half way between the charges, where its magnitude is 


ttT 


(4.22) 


Away from this line the exponent in (4.19) enforces an exponential decay in the 
magnitude of the held. Indeed we hnd 


(x-a;^ (tq))^ 
2G (to, To) 


|x - a| |x - b| (|x - a| + |x - b|) 

/ , , , , |x — a| + |x — b| 

-2 (x - a) ■ (x - b) ^ 


|x-b| 


X — a| |x — b 


(4.23) 


and the direction of the held depends on the spatial point through the hrst term 
in (4.17). It is possible to use (4.19) to determine the contribution at an arbitrary 
spatial point but the result is less illuminating than a visual representation of the 
held lines and the held magnitude. It is most useful to plot the streamlines, tangent 
to the held at each spatial point. Fig. 2 and Fig 3 show the held strength and 
direction on the z = 0 plane of a pair of oppositely charged particles placed at 
positions a = (3, 0,0) and b = (—3, 0, 0) for two values of the parameter T. We have 
imposed a sharp cut-oh about the positions of the charges to avoid the divergence 
encountered there. 
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(a) The magnitude of the electric field (b) The stream lines of the electric field 


Figure 2: Field magnitude and field streamlines in the low T limit - in this plot 
T= ^|b-a|^ 

oO ' ' 



(a) The magnitude of the electric field (b) The stream lines of the electric field 

Figure 3: Field magnitude and field streamlines in the low T limit ~ in this plot 
T = |b — a|^. The field decays exponentially as described in the text which 
gives rise to the white-space in the plot, in which the field is negligibly small. 

The form of these field configurations suggests that the low T limit of this theory 
gives some sort of confining field, albeit not one with a potential linear in the separa¬ 
tion of the charges. In this way T interpolates between the classical field associated 
to a pair of charges and a regime in which field lines are concentrated about the line 
joining the charges. The emergence of confining behaviour has briefly been consid¬ 
ered for action at a distance theories before [42, 43] and we add to the discussion 

®This expression is consistent with the behaviour of the system under exchange of a and b 
because this is equivalent to sending r —^ 1 — r. 


25 









with our determination of the held lines with this new model. The three dimensional 
plot in Fig. 4 highlights how the low T lines of hux are compressed into a thin tube. 
This completes our discussion of the correction to the bosonic theory at hnite T. 



Figure 4: Conhning held lines in the low T limit - here T 


1 

500 


b- 


|2 

a| . 


4.2 Corrections in the fermionic case 


In this case we shall consider four dimensional space-time and will briehy state the 
result of analysing the order ^ corrections to (3.30). The general approach is identical 
to that of the bosonic theory. We have already determined the momentum space 
correction to the hrst and last terms in (3.27) in the previous section. For the 
contribution corresponding to the end of the curve at the fith component is given 
by 

2q 
T 




k^k’' 

k^ 


{h-a) 


Ak-a 


(4.24) 


We also hnd a subleading contribution from the middle term in (3.27) which takes 
the form 

^ (i“ “ (4.25) 

SO now we must determine the four dimensional integral 


(Ak 

(27r)^ k^ 


ik-a—x 


(4.26) 


Symmetry dictates it must be equal to + B {a — xY {a — xY |a — x\~‘^ where 
the constants A and B have dimension [length]”^. Proceeding as before we contract 
with and {a — xY A — xY to produce two equations. The hrst gives 


AA + B 


(Ak ^ _i 

(27r)"^ k"^ \x — a 


(4.27) 
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and the second leads to 


la - (A + B) = j ^ 


(4.28) 


This integral can be carried ont by hrst doing the integral with respect to where 
donble poles are enconntered at ko = |k|. The remaining integral over the three 

dimensional vector k can be done by choosing the z-axis to align with x — a as above. 
Then we are left to determine 
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(4.29) 


Integrating over 9 leaves only 
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(4.30) 


Carrying ont the fc-integral several times by parts yields \x — a'f {A + B) = ^ so 
(4.26) evalnates to 


Stt^ 


X — a\ 


5^’^ -2 


(a — x)^ {a — x) 


X — a\ 


(4.31) 


So in position space (4.24) becomes 


-Q 


\x — a| 


{h — aY + 2{x — a) 


^ (x - a) ■ (6 - a) 


\x — a\ 


(4.32) 


for which it can be verified that the divergence vanishes, and (4.25) becomes 


-q 


Ati^T \x — a| 


(5^^ - 


{b — aY — 2 (x — aY 


(x — a) ■ (x — a) 


X — a 


(4.33) 


which is also divergence free. The first order correction is fonnd by snbtracting these 
and inclnding the contribntion from the other end of the cnrve. We hnd at order P 
the expectation valne of the insertion evalnates to 
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The form of this correction has a similar functional form to that of the bosonic 
particle discussed above. 

We may also ask about the low T expansion to determine the behaviour of the 
system in this regime. It is again useful to carry out the integral over k hrst to arrive 
at 


(F {x))^ = 


-q dri 

{27rfJo (TG{Tun)y 


2G (ri,ri) 


[X 


- UJcY - 




{x - UcT 
2G (ri,ri) 


( 5 ^"^ - 7 ^ 7 "^) 


exp 


{x-UJcf ^ 
2TG 


(4.35) 


In the limit that T vanishes we have a representation of the four dimensional function 
supported on the straight line joining the charges in Minkowski space. Applying 
Laplace’s approximation allows us to determine the leading order behaviour which 
we illustrate for a point x on the line from a to b. In this case we clearly hnd the 
maximum at the point where Uc (tq) = x (tq remains unchanged from (4.20)) so the 
hrst and last terms vanish. Explicit application of (4.19) results in 

q h — a 

(^27rTG(ro,ro))" 

= —^- z{b-a) (4.36) 

(27rT)2 |x-a |2 \x-h\~^ 


which achieves a minimum at the midpoint of the line of magnitude 

q 


(4.37) 


For sufficiently small T, this will be greater than the corresponding minimum asso¬ 
ciated to the bosonic theory. 

This section contained some analytic results for bosonic and fermionic point 
particles beyond the leading order behaviour. It is especially interesting to note 
the small T limit of the system which demonstrates a localisation of the held about 
the classical path between the charges. In the following section we return to the 
dehnition of the interaction between particles and use it to construct a full quantum 
theory. 


5. Contact interactions between particles 

In this section we extend our work to describe a set of particles which interact when 
their worldlines intersect. Accordingly we work in a D = 4 dimensional space- 
time. This section carries out the analogous analysis to that in [2] where a collection 












of strings interacting upon contact was considered. In this section we limit our 
discussion to bosonic particles. To describe the dynamics of a collection of interacting 
particles we augment the free action, So, of each point particle with a non-local 
contact interaction as follows: 


‘S'tot — So [ci, uji] + ^ 





dTidTj Ui (n) ■ Uj {tj) (Ui (Tj) - Uj {tj)) (5.1) 


There are two ways to motivate this interaction. In correspondence with the intro¬ 
duction of [2] one can substitute the expression for E', the insertion used in Section 
2, into the standard action of Maxwell electromagnetism for the energy of a static 
electric field: 


Sf 


= ^ y d^x E ■ E = ^ y y dTidT 2 Co (ri) ■ Cc {T 2 ) 5^ {uj (ri) 


lj{t2)) (5.2) 


which produces a self-interaction along the worldline in question (we have integrated 
over the spatial point x). It is natural to then extend this to encompass interactions 
between distinct particles as we have done in (5.1). An alternative is to follow the 
worldline approach® and consider two particles which interact via the exchange of 
a single massive gauge boson between the points and The propagator for 
this process may be taken as { 001 , 002 ) oc y J d'^p^—ASA^, -which (as has been 
discussed by Feynman [44] - see also [7]) should be inserted into the path integral 
and integrated with respect to the points of emission and absorption along each 
worldline. To take into account all possible interactions of this type, one is thus 
required to include a term 


JJ duidu2A^^{ui,U2) (5.3) 

in the worldline theory. To achieve (5.1) from here it suffices to parameterise the 

paths by variables ti and T 2 and take the limit that the mediating boson has large 
2 

mass so that ^ 0. This provides the appropriate (5-function contact term, since 

the range of the interaction will be heavily supressed in this limit, and the resulting 
expression can be identihed with (5.1) by setting = gm"^. 

The form of (5.1) may appear unusual but it is straightforward to verify that 
it satishes the consistency criteria described by Kalb and Ramond for interactions 
of this type [18, 17]. We also recall from the discussion in the introduction that 
the (5-function ensures that the interaction is local in space-time. We must specify 
fixed boundary conditions for each particle, which we denote by ooi (0) = af and 
(Uj (1) = The coupling strength g determines the relative strength of interactions 
between the worldlines (we have taken the couplings between all pairs of particles to 

am grateful to the reviewer of this article for their insight in pointing this out. 
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be the same, rather than introducing separate coupling constants gi as presented in 
the introduction) and its sign determines whether the interaction is repulsive {g > 0) 
or attractive < 0). 

We shall consider the partition function of the theory described by Stot and 
determine its physical content as well as investigating whether there are divergences 
which need regularising. In order to build upon our earlier results which allowed us 
to relate the worldline theory to electrostatics we will retain the constraint on the 
path length of each worldline which was used in earlier sections, so that for each 
particle we impose f dTiei{Ti) = T. Since there is no reason to keep hold of this 
arbitrary dimensionful parameter we will follow the choice made in sections 2 and 
3 by considering the large T limit. We expect this limit to be the correct one for 
(5.1) to provide an extension of electrostatics and to ensure the interaction between 
particle worldlines is not unduly restricted to a localised region of space-time. In this 
theory we do not anticipate divergences corresponding to the coincidence of operators 
due to the well behaved nature of Green functions on one dimensional domains but 
we do stand to encounter unwanted divergences in taking the T —)■ cx) limit. We will 
also hnd that the dehnition of the contact interaction will require slight rehnement 
corresponding to the rather trivial vanishing of the argument of the 5-function when 
T = t' and the worldlines Ui and uj are the same. 

Before proceeding we pause to point out how this action differs from the quantum 
mechanics which arises in both the action at a distance theories and the worldline 
formalism of a scalar quantum held. To compare with the action at a distance 
theories discussed in the introduction we may dehne a functional of the worldlines 
in order to recast the action above as if the particles interacted with a current. Let 



(5.4) 


which we know from section 2 satishes (x) = 6^ (x — Qi) — 5^ {x — bi) so is 

sourced only at the points corresponding to the end of the worldlines (this is necessary 
in order that they be held in place). Then the interaction between the particles can 
be written as 



(5.5) 


and the equations of motion for the particles take the form (for illustration we gauge 


hx the einbein to e = 1 which imposes the hrst class energy-momentum constraint 




(5.6) 
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which are well-known for the minimal vector coupling between worldlines and a 
background current^®. We must recall, however, that this background current is 
not an independent degree of freedom because it depends on the particle worldlines 
(which must satisfy the equations of motion) and is also not a conserved quantity. 

We also comment briefly on the relationship to the worldline formulation of 
quantum held theory to further justify its study. Although we have discussed how 
the interaction can be thought of as arising due to the exchange of very heavy gauge 
bosons between the worldlines, it can also be compared to the transition amplitude 
of the worldlines in some potential. If we consider a massless scalar held, 0, with an 
interaction potential U (0) (for simplicity we suppose that 0 does not couple to any 
gauge held) then the one-loop ehective action, T [0], is dehned by integrating over 
the quadratic quantum huctuations of the matter degrees of freedom [45]: 

(5.7, 

This is then rewritten in the worldline formalism by expressing the logarithm as 
an integral using the Schwinger proper time trick [46] and the functional trace as 
a quantum mechanical transition amplitude for a point particle to traverse a closed 
loop: 

Trln oc ^ (^| exp [-T [-d'^ + rri^ + U" (0 (cu (r)))]) |a;) 

(5.8) 

= dr^^+m'^T + U" 

The ehect of the self-interaction potential is to modify what would otherwise be the 
quantum theory of a free particle via an additional term U” (0 {u (r))) in the action. 
In the case that m = 0 and there is no held potential the particle action reduces 
to So which we used for functional averages in section 2. There we carried out the 
integral over T by making use of the constraint on path lengths. 

The calculation of the ehective action must now proceed perturbatively (for 
example, by expanding 0 about the centre of mass of the worldline, which is the 
zero mode of the kinetic operator [8]). However, if for the sake of illustration we 
wish to calculate the 2A^-point one-loop scattering amplitude, we may functionally 
diherentiate Wtimes with respect to 0 and then expand the held as a sum of N 

^°For completeness we point out that the equation for the “field” can be imposed by intro¬ 
ducing a trivial space-time action | f d'^xJ^J^ whose variation, in conjunction with that of (5.5), 
ensures the correct on-shell assignment in (5.4). This is slightly different from the form of the 
theory discussed in the introduction where the interaction is truly non-local and is carried by a 
Green function with may be supported away from the null light cone. In that case the equations of 
motion and the space-time action would be familiar from minimal electromagnetism. 
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plane waves: 


N 

= (5-9) 

i=l 

Now, the physical picture of the interaction in (5.1) is of a pair of worldlines meeting 
at a point, where they interact before independently separating. This reminds us of 
the Feynmann diagram in 0“^ theory so let us consider U (0) = |t 0^, which leads to 
an additional worldline interaction |0 {uj (r))^. Carrying out the Fourier integration 
over momentum the one-loop 2iV-point amplitude (0 (xi)... 0 {x 2 n)) can be written 
in position space as a sum over permutations of pairings of the form 


N 




dTT 




n / (f) - X2i) {u {Ti) - X2i-l) e 


^2 

0 dr ^ 


2=1 


( 6 . 10 ) 

The interpretation of this is that the ^-functions force the joining of pairs of points X 2 i 
and X 2 i-i to a point on the virtual worldline oj (tj), reproducing a quartic interaction 
vertex. The interaction that we propose in (5.1) is slightly different since we couple 
pairs of worldlines to one another directly without the need for a virtual loop to 
mediate the process but the interaction when pairs of lines intersect has the same 
form. The one-loop example considered here can be interpreted in terms of worldline 
interactions by imagining a particle travelling between two space-time points whose 
worldline intersects a second, closed worldline at a single point only. 

Our integrand, however, does not consist only of the minimal scalar measures 
J dxi but carries with it directional information^^ through the factors of So we 
imagine a four point contact interaction between pairs of worldlines which depends 
also on the tangent vectors to the worldlines at the point of contact. We get an idea 
of what this would mean for the corresponding held theory by dimensional analysis 
of the coupling strength. In Zi)-dimensions, [g] = length^”^, whilst a scalar held 
has [0] = length If we introduce a coupling constant A then the held-theory 
interaction ^ f d^cp 0 ^ d^x is dimensionless if the dimensions of A are 

the same as g. This suggests that the worldline theory proposed here is related 
to a held theory with quartic coupling containing two derivatives (we must recall, 
however, that the constraint on path lengths being imposed means that we are not 
reproducing the entire content of the held theory, which would require us to integrate 
over T rather than consider only the contribution as T —)■ cx)). With this heuristic 
picture of how the contact interaction can be interpreted we now return to determine 
the partition function of the modihed worldline theory. 

We shall carry out the calculation as a perturbative expansion in g. Expressing 


wish to thank Paul Mansfield for helpful discussions on the physical interpretation of the 
interaction described in this paragraph. 
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the (5-function in its Fourier representation introduces vertex operators 


-^Vk {Ti) ■ V_k {Tj) (5.11) 
(27r) 

At first order in g the correction to the partition function of the non-interacting 
theory takes the form 

1 SidTi -T^ g-So[aM,l'j jj dujj- dUi S‘‘ (uij - Ut) 

(5.12) 

There are two contributions to this sum. When the worldlines are distinct (5.12) can 
be factorised to make use of the result in section 2: 
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3 So [ej ,ujj 


duj^ 



du^d'^ {cuj - Uk) 
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(5.13) 




GjdTj- 


^-So[ej,ujj] 




j \ \^j I 



(5.14) 


where we have taken the large T limit and discarded contributions of order ^ arising 
from the expectation value in the hrst line. The integral over ooj at the far right of 
the bottom line produces a boundary contribution and we recall that the boundary 
conditions ensure this is a constant throughout the functional integral over ojj. So 
the functional integral is rendered trivial and the result is 


—y 

4^2 

j¥=k 





(5.15) 


The interpretation of this is a pairwise interaction between the hxed ends of the open 
worldlines of the theory, whose form we recognise as the space-time Green function 
(d^) This is the Green function of the scalar held theory whose worldline reformu¬ 
lation leads to the quantum theory described by Sq which we have discussed above. 
The effect of the contact interaction is to couple the end points of the worldlines 
to one another via the propagators in (5.15). This retrospectively suggests that the 
high T limit was a good choice, since it leads to a familiar interaction with the long¬ 
distance properties that one would hope for in quantum held theory, in much the 
same way that the same limits taken in sections 2 and 3 yielded familiar electrostatics 
helds. 
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This can be contrasted with the result we found in the case of interacting strings 


in [1] and [2]. We have already explained how the endpoints of the strings were 


fixed to lie on some specihed worldlines. The effect of the contact interaction on the 
worldsheet came primarily from the region close to its boundary, B, and provided 
an interaction between the worldlines which took the form 



(5.16) 


which now contains the space-time propagator of a vector field integrated over all 
points on the boundary of each string. In fact the Fourier transform of the Green 
function took the form (k) = — k^k^/, which is the propagator for 

an Abelian gauge boson, A in the gauge d ■ A = In this way we related the 
tensionless limit of the interacting string theory to quantum electrodynamics, where 
the contact interaction served to produce the interactions between the worldlines 
and the gauge field. The calculation of the partition function of the point particle 
theory is related to (5.16) as analogous to the contribution arising from two fixed 
points in space-time. If these are taken to be points oJa and Uh on two curves Ba 
and Bb then the string theory result arises by then integrating over all such points 
on the curves (see also the end of section 2 where we discussed the transition from a 
spatial average to averaging over curves in space-time). In the string theory version 
this integral is generated by the extra structure in the worldsheet contact interaction 
which produces the integration measure and automatically counts the contributions 
from every point on the string boundary. 

Returning now to the point particle theory we must address the case j = k 
separately as then the worldlines being integrated over are the same. In the string 
theory case [1, 2] we found that two contributions make up this interaction - a 
renormalisation of the free string action and cosmological constant and a contribution 
corresponding to self-intersection of the string. For point particles, however, no such 
self-intersection ought to be present since the vanishing of 



(5.17) 


is only at Ti = T 2 - Naively this would provide (0) x Length (ca), suggestive of 
a renormalisation of the non-interacting part of the action. A second major differ¬ 
ence between the theory of point particles we consider here and the string theory 
equivalent is that the current case does not require us to worry about encountering a 
conformal anomaly. In particular there are no short distance divergences associated 
to the Green function which would require regularisation. To make all this more 
precise we must consider 



T 


(5.18) 


34 






where we have again used the Fourier representation of the four-dimensional 6- 
function. The insertions are once more easily generated via the introduction of 
sources and the integral over u provides the following generalisation of (2.12) 


9 [ f f, , S 6 

2] (27r)^yy {ri) {t 2 ) 


exp 



■ J {Ti,T) 


kj G {Ti, Tj) — i dr k, 

9 i=i 

exp (- j drj (t) ■ (t) + ^ ■ Uc (n) j 


j=0 

(5.19) 


where ki = k = —k 2 - It is trivial to carry out the functional differentiation and it 
proves useful to dehne 


'I' ('^1, T 2 ) = G (ti, Ti) -|- G {t2, T 2 ) — 2G (ti, T 2 ) (5.20) 


in order to express the answer as 


I / dn dT 2 exp (ri,r 2 ) 




+ -k^u^^dl^ (ri,r2) - - - 


1^ 


dtdt'^ (ri,r2) 


^ik-(u)c{T{)-Ldc{T2)) 


(5.21) 


which must then be integrated over k. In the above equation we have integrated by 
parts to produce the transverse projector for the hnal term in square brackets. In the 
worldline formalism it is often advantageous to follow an alternative integration by 
parts procedure advocated by Bern and Kosower [47, 48] which removes all second 
derivatives of Green functions^^. The interpretation of this technique has its roots 
in the pinching of Feynman diagrams present in the underlying held theory, but as 
we do not necessarily have such a model behind our work it is unnecessary for us to 
adhere to it. Now 

{n,T 2 ) = T (|ri - r 2 | - (ri - T 2 f) (5.22) 

actually coincides with twice the Green function with periodic boundary conditions 
[38]. As such it satishes (b, ^ 2 ) = d {ti — T 2 ) — 1 and is a function of ri — T 2 

^^The worldline approach was originally proposed by Strassler as a means of arriving at earlier 
results relating field theory amplitudes to the infinite tension limit of string theory [49, 50], so that 
many calculational tricks used by string theorists are borrowed for worldline calculations. 
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only. Noting also that Uc (ti) — Uc {T 2 ) = {b — a) (ri — T 2 ) we learn that the entire 
expression (5.21) is in fact a function of the separation ti — T 2 - We use this to hx 
the zero by setting T 2 = 0 and multiplying by dT 2 = 1. Furthermore d' (ri, 0) = 
Tti (1 — Ti) is just the coincident Green function we met in section 2 so we must 
determine 


g J dn exp G (n, ri) 


in, Ti) + ^-k'^u^dlG {n, n) 


ri‘ 




k^k’' 

1^ 


dtdt^ in, 0) 


Jk-{b-a)Ti 


(5.23) 


Anticipating that we will eventually take the limit T —)■ 00 we need only consider 
an expansion in powers of The leading order contributions again come only from 
the ends of the interval where the coincident Green function vanishes so it suffices to 
expand about ri = 0 and ri = 1. We shall consider each term in (5.23) separately; 
at lowest order in ^ the hrst takes the form 


= % - 1 ] ^ + O (^) (6.24) 


Upon integrating over k the latter term in square brackets vanishes whilst the former 
provides (27r | & — a|)~^. The contribution we seek is the trace of this - we note that 
since the derivative of the classical solution to the equations of motion is a constant 
we may express the above expression as 


9 

\b — a f 



(5.25) 


We interpret this as providing a renormalisation of the free action and note that it 
is suppressed in the large T limit. 

The second and third terms in (5.23) involve a derivative of the Green function 
and as such provide contributions that are independent of T. Taking the trace we 
require 


gik ■ die 


J dn G (nn) e- 5 ^^G(n.n)gifc-(&-a)n 


tkjJ^ r ik.{b-a) 

^ 2fc2 



Garrying out the integral over k provides 


Cjc- {b- a) 
\b — a\^ 


+ 0 



(5.26) 


(5.27) 
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Recognising that u = b — a we may also cast this into the form 


9 

\b — a 



+ 0 



(5.28) 


providing a finite renormalisation of the free action. 

Finally we take the last term of (5.23) and consider its trace. We nse the defining 
eqnation of the Green fnnction to write its contribntion as 


-g / dTi T (b (ti) - 1) e 2 


-^k^G(Ti,Ti)^ik-(b-a)Ti 


(5.29) 


Both terms are independent of the field Uc which snggests that we shonld interpret 
them as renormalisations of a cosmological constant term in the action dr T. The 
(5-fnnction gives ns 


-9 drT 


(5.30) 


Jo 

exactly as reqnired, albeit formally divergent when we take T to inhnity^^, whereas 
the second evalnates to 





which npon integrating over k becomes 


_9 _^ 

27r2 \b- a 


+ 0 



(5.33) 


(5.34) 


This snrprisingly takes the same form as (5.15) and snggests that there is after all a 
self-interaction present in the theory, sensitive only to the bonndary of the worldline. 

To lend this more weight we conld approach the calcnlation in a complementary 
fashion. We note that T (ri,r 2 ) vanishes only at T 2 = ti so we conld arrange onr 
calculation by instead expanding about this point. When ti is not close to the 
boundary (measured with respect to ^) integrating T 2 about ti corresponds to 

^^The integral over k also provides an infinite volume factor multiplying this result. We may tidy 
this up by carrying out the fc—integral before integrating over ri. Doing so turns (5.29) into 

J ( 2G(r,.n) ) ' 

We are concerned only with the contribution arising from the (5-function so we need the value of 
the exponent as ri —>■ O’*". An easy calculation shows that the exponent vanishes leaving only 

-9 (0 " (G (ri, ri))"T (5.32) 

which diverges. 
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integrating over their relative separation and the leading order contribntion arises by 
setting the T-independent exponent exp {ik ■ {h — a) (ri — T 2 )) eqnal to nnity. This 
gives the terms above which are absent of an exponent. When ti is close to the 
bonndary we must take care because T 2 is restricted to lie in [0,1]. So for example 
when Ti ~ 0 we must integrate T 2 a small distance from this boundary along the 
line but must also consider the contribution when T 2 is integrated from the opposite 
boundary along the line. Indeed for the latter case we consider, following the notation 
of section 2 , 

1 / dn f dT2 T {5 (n - T 2 ) - 1) ^5 35 ^ 

2 Jo Jl-h 

The h-function is not supported for this configuration of variables, besides it is the 
second term with which we are concerned. At leading order in ^ we evaluate the 
trailing exponent at the point ri = 0, r 2 = 1. A change of variables r = 1 — r 2 yields 
T (ti, 1 — r) = T ((ti + r) — (ti + r)^) and the integral becomes 


I'h /‘2h 

-- dn du 
2 Jo Jo 

where we have set u = ti + t. At leading order this provides 

9 —ik-(b—a) 


(5.36) 


(5.37) 


which is to be compared with (5.33). The sign of the exponent is not important 
and the factor of two is found by including the configuration where the positions of 
Ti and T 2 in the above calculation are swapped. So the physical information arises 
when the two points are close to opposite boundaries, whereas the renormalisations 
appear from the region where the two points coincide. This is in fitting with the 
naive analysis of the form of the contact interaction when the worldlines are the 
same. 

The self-interaction (5.34) also has the appropriate factor of two difference for 
it to be subsumed into the sum (5.15) so that at first order in the expansion of the 
interaction we hnd 


9 

3,k 





(5.38) 


We write to denote that when j = k we discard the hrst and last terms in 
the summand where the separation vanishes. This concludes the analysis of the 
hrst order effect of the contact interaction present in this theory. In the following 
section we turn to consider higher order interactions and show that the result above 
is not spoilt and that this propagator interaction between the end of the worldlines 
exponentiates. 
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5.1 Expansion to arbitrary order 


We may consider expanding the part of the action corresponding to the inter-particle 
interaction to arbitrary order in the coupling strength. When the interaction is 
between distinct worldlines we may repeat the analysis at hrst order to simplify the 
calculation to a form familiar from section 2. Here we address the problem of having 
a hxed even number, N, of points to be integrated over the same worldline which 

N 

could potentially spoil the result at order or higher. The reason for this is that 
we must be cautious in taking the limit T —>■ oo, in case the clustering of points on 
the worldline produces dependence on T which would cause the result to diverge. So 
we consider 


K (^i) K ■ ■ ■ K 


(5.39) 


where we understand that each point must be integrated from r* = 0 to 1. 

Wick’s theorem produces a common factor to all of the contractions that could 


be formed out of the above product of helds which takes the form 



(5.40) 


exp 


We are interested in taking the high T limit, whereby the leading contribution will 
be from the regions of integration where the exponent vanishes. There are no short 
distance divergences which require regularisation so we may consider the points to be 
arbitrarily close to one another. This leads us to consider two cases. The exponent 
can be made to vanish by ensuring that each individual Green function G{Ti,Tj) 
vanishes, which requires us to localise each of the points close to either end of the 
domain. Alternatively we might consider bringing a collection of points coincident, 
which is instructive in order to compare with the results of the previous subsection. 
As no regularisation is required we keep the discussion of this case brief. It is useful 
to split (5.40) as 



(5.41) 


which makes it clear that for an arbitrary conhguration of the points away from the 
boundary the hrst factor damps the integrand. Now suppose that some number, n, 
of these points are brought to the same point ti. Then these points contribute 



(5.42) 


exp 


where the ... represents terms which depend on the relative separation between each 
point and ri. Recall that G{ti,ti) = TG{ti,ti) vanishes only at the boundary of 
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the domain. Consider how this exponent behaves under the integral with respect to 
fci, say. In the large T limit a Laplace approximation implies its effect is to provide^'^ 


(2rG(r„Ti))" 


(5.43) 


Since in one dimension the coincident Green function is finite we see that all contri¬ 
butions from such a configuration are suppressed by a factor of T“^. The size of the 
integrals over the relative separation r* — Ti can be examined by power counting in 
The largest contribution arises when all fields take part in a contraction and can be 
arranged into a series of | second order derivatives. A simple calculation shows that 
this contribution is of order Tt. But after integrating the other n — 1 points about ti 
the resulting expression is of order T. In combination with (5.43) this contribution 
can be seen to be subleading in 4^ so no unwanted divergences are encountered. 

Returning to the case that all of the points are on the boundary, each of the Green 
functions in (5.40) vanishes. So we expect a contribution to the expectation value 
from the region of integration where each point is close to the boundary. We have 
learnt in the previous subsection that the O {!) term arises when we contract the fields 
inside each vertex operator amongst themselves. The exponents exp (ifc ■ cuc (t*)) 
can be approximated at leading order by replacing them with their values at the 
appropriate boundary - we denote this by exp {ik ■ ujeJ- So we consider a term with 
r such contractions: 



The leading order contribution clearly requires r = N so that all fields are contracted, 
whereby the above expression reduces to 


n 



Akj-a 






(5.45) 


We now impose pairwise /cj+i = —ki and contract /Xj and /ij+i to produce the expec- 

N 

tation value of the contact interaction at order : 


N 

-2^i 


2 ik-{b—a) 


n 

2=1 


kf 


(5.46) 


^^This should be considered in light of the results at first order where we had ki = k = —k 2 
where the exponent vanished throughout the domain when ri = r 2 . We exclude such cases here 
because we understand that they lead to a simple renormalisation of the action. 
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As our analysis has been for an arbitrary number of points the result at this order 
is not spoilt by considering a higher order term in the expansion of the contact 
interaction. The above equation is simply the Fourier space version of the order 

contribution to the exponential of the sum of boundary interactions (5.38), from 
which we conclude that this interaction between worldlines exponentiates absent of 
short-distance divergences. 

6. Discussion 

In this paper we have considered contact interactions in the context of theories of 
point particles. We have demonstrated a novel way of generating the static electric 
held for a pair of point particles by considering huctuating worldlines whose end¬ 
points are hxed to the positions of the charges. The functional approach we used 
for calculation allowed us to generalise this construction to include spin 1/2 particles 
and we then used the formalism to construct a quantum theory describing a set of 
point particles interacting upon contact. The result for fermionic particles is inter¬ 
esting because it suggests an unusual electric force acting on Dirac spinors due to the 
charged particles. Both results depend somewhat on the choice of boundary condi¬ 
tions imposed on the worldline helds at either end of the paths and we also found it 
necessary to introduce constraints on the worldline metric and its super-partner to 
reconstruct the classical held of static charged particles. 

We demonstrated that the worldline contact interaction provides an unconven¬ 
tional method of generating the static dipole electric held due to equal and oppositely 
charged particles, and also modihed the boundary conditions on the curves to pro¬ 
duce the held of a static point charge. This result required us to take the intrinsic 
lengths of the curves, T, to be very large but we also investigated the corrections to 
the helds for hnite values of T. The inverse square law for electrostatic interactions 
has been verihed to a very high precision which provides a bound on the size of the 
deviations from this behaviour. This in turn allows us to hnd a minimum bound for 
T. An alternative method of estimating the size of T follows from the 1/T correc¬ 
tions to the retarded Green function given in equation (2.30) which are discussed in 
section 4. This would imply corrections to the speed of propagation of light, c, which 
is also known to extremely high precision. 

Should future experiment further rehne these two quantities then we may in¬ 
crease our lower bound on T, or we may test unexpected deviations against the 
results provided by the model in this paper. Either way, from a theoretical stand¬ 
point varying the parameter T has an important effect on the form of the held lines, 
especially when we consider the low T behaviour where the novel behaviour of con- 
hning held lines was discovered. It is important to better understand this limiting 
case because it may help to develop another approach to this well known aspect of 
non-Abelian gauge held theories. The worldline techniques we have used here may 
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prove valuable in approaching this problem because of the computational simplicity 
present with this approach. 

Beyond the hrst order calculations we presented in the hrst few sections of this 
paper we also determined the partition function for a set of particles which interact 
when their worldlines intersect. We showed how the contact interaction produces 
propagator couplings between the hxed endpoints of the worldlines, which is not 
unexpected considering that the effect of the interaction is to modify the underlying 
scalar held theory to include a quartic self-coupling. An important avenue for further 
investigation would be to work out the exact relationship between the quartic held 
theory and the interacting worldlines, since it may then prove possible to rephrase 
old questions in held theory in terms of simple quantum mechanics on the worldlines, 
adding a calculational tool which readily inherits the efficiencies present within the 
worldline formalism of quantum held theory. 

The worldline formalism highlights an intimate connection between held theory 
and the hrst quantised particles we have considered here. Our approach dihered in 
that rather than coupling the theory living on the worldline to a background gauge 
held (as occurs when integrating out matter helds with the worldline approach) 
we coupled the worldlines to one another directly using the 5—function interaction. 
We also introduced fundamental constraints on the worldlines’ lengths (rather than 
integrating over T) and worked only with open worldlines. However, the functional 
methods that we used are not limited to the calculations we have undertaken in this 
instance and can be directly applied to other situations that arise in the worldline 
approach. 

It would be interesting to consider a similar contact interaction dehned on closed 
worldlines since this case is encountered in the determination of one-loop ehective 
actions and multi-loop scattering amplitudes in the worldline approach. There are 
substantial differences between the closed worldlines and the open paths that we have 
considered in this paper, in part due to the difference in their topology. In particular 
the Green function is related to T(ri, T 2 ) of (5.22), the boundary conditions requiring 
it to be periodic rather than vanishing at the ends of the interval. In this case, 
the Green function’s coincident limit is constant along the worldline (momentum 
conservation implies that, in flat space, the numerical value of the coincident limit 
is not important). In this article we relied on the functional form of the Green 
function to ensure that in the high T limit our integrands were suppressed moving 
away from the endpoints of the worldline which would no longer hold true if the 
particle’s path were closed. It would thus be necessary to improve the analytical 
calculation of the integrals instead of extracting the leading order behaviour as we 
have in previous sections. At higher loop orders the worldlines have extra insertions 
of propagators and so we would like to generalise this work to the case of an arbitrary 
one dimensional topology (the Green functions for such worldlines have been studied 
extensively in [38, 51]). This would signify useful progress in the hrst quantised 
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formalism. 

We have also discussed how the calculations presented here for worldlines with 
hxed endpoints are related to an analogous theory of interacting strings, where the 
worldline formalism allowed us to relate the string theory to products of Wilson 
lines in spinor quantum electrodynamics. As in this paper, it was also necessary to 
take a particular limit in order to secure that result, namely the limit of vanishing 
string tension. This has the effect of making the strings macroscopically large in 
comparison to the worldlines (which are their boundaries), just as the T —)■ cx) limit 
discussed here makes the worldlines’ lengths very large. In contrast to this paper 
the tensionless limit of the string theory was necessary to suppress divergences which 
would otherwise spoil the conformal invariance of the classical theory, which presents 
a fundamental obstruction to investigating hnite tension corrections. In this sense, 
the worldline theory has an advantage over the interacting string theory, in that we 
may vary T without encountering any quantum anomalies. 
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Appendix A - Evaluating the free particle determinants 


Here we determine the normalisation constants used in section 2 using (^-function 
regularisation. With our gauge choice (e = T) we make a change of variables t = Tt 
and dehne 


Z' = 


faj(T)=b 


rT j 

^(jj e -^0 2 


J a>(0)=a 

which can be interpreted as the matrix element (b 




( 6 . 1 ) 


a ) with Hamiltonian l-L = 


In the Schodinger representation this becomes the position space representation 


of the heat kernel: ( b 


_ rj-t V 

6 2 


ay. This is a well know expression but we hnd it using 
functional methods in keeping with the spirit of the remainder of this article. 

Generalising to arbitrary dimension D, the integral over uj gives a functional 
determinant and a boundary term: 


= ldet(T 


(b-ar 

e 2 t 


( 6 . 2 ) 


The determinant of an operator, O, can be dehned using the (^-function [52] as 


det (O) =exp ( 


2=0 


Co (z) 


(6.3) 


n=l 


where the A„ are the eigenvalues of O. The formal product of eigenvalues is then 
regularised via the analytic continuation of the (^-function. With Dirichlet boundary 
conditions the eigenvalues of the operator in Z are A„ = (^) so 

(6.4) 

which has derivative C' (0) = 21n (y)C (0) + (0). So with this regularisation we 

arrive at 


giving 


, D (b-a)2 

Z' = (27rT)" 2 e— 


( 6 . 6 ) 


If instead we consider mixed boundary conditions uj (0) = 0, ch (1) = 0 then the 
eigenfunctions are 

^ (6.7) 


sm 


where the eigenvalues are now given by A„ = ^ 
obtain 


For the (^-function we now 

( 6 . 8 ) 
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where we’ve introduced the Hurwitz zeta function ({s,q). The derivative with re¬ 
spect to 2 ; is then (' (0) = 21n (^)C (O, + 2^' (O, Now ( (O, |) = 0 so the hrst 

term vanishes along with the T dependence, leaving 


det 



= 2 


(6.9) 


for the case of mixed boundary conditions. The change in boundary conditions also 
alters the boundary contributions from the classical action and we hnd that Z' is a 
constant independent of T and a. 

For fermionic helds the kinetic term is hrst order in derivatives and we impose 
anti-periodic boundary conditions in the case of closed paths which represent traces. 
We may of course calculate 

Zf= [ e-^o V’-Wr ^ 

i^(r)=-i/.(o) V \dtJJ 


using standard techniques from quantum mechanics but for completeness we continue 
to apply functional methods. We need the eigenvalues A„ = and form the 

determinant via 


n 


(2n -|- 1) Tii 
T 


( 6 . 11 ) 


n=—oo 


The (^-function for this operator is thus 



— 2 


c 



( 6 . 12 ) 


which has derivative C (0) = In (^)C (O, + C (O) |) which evaluates to — In \/2 so 


ZiT = 2f 


(6.13) 


which we note does not depend on T ~ indeed the change of variables to r G [0,T] 
does not change the form of the kinetic term. 

The other result we need is the normalisation constant for the open path action 
which includes the term dr ^ui ■ -0. We use (3.16) to calculate 

J dxoHxo-^)^i^ ^ J dxoHxo - ^ {e^ . (6.14) 


Integrating over xo picks out 


1 + 


Xb 


7 


2V2T 


which is to be multiplied by Zp calculated above. In the main text we impose 
which vastly simplihes the remaining calculations because it decouples u and 


(6.15) 

S = 0 
xp. 
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